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O ■ Abstract 

o , 

CsJ \ The purpose of this work is the computation of the cosmic Type la supernova rates, namely the frequency of Type 

^ ■ la supernovae per unit time in a unitary volume of the Universe. Our main goal in this work is to predict the 

Type la supernova rates at very high redshifts and to check whether it is possible to select the best delay time 
distribution model, on the basis of the available observations of Type la supernovae. We compute the cosmic Type 
la supernova rates in different scenarios for galaxy formation and predict the expected number of explosions at 
high redshift (z > 2). Moreover, we adopt various progenitor models in order to compute the Type la supernova 
rate in typical elliptical galaxies of initial luminous masses of IO^^Mq, 1O"M0 and IO'^Mq, and compute the total 
amount of iron produced by Type la supernovae in each case. In this analysis we assume that Type la supernovae 
Q,! ^s caused by thermonuclear explosions of C-O white dwarfs in binary systems and we consider the most popular 

frameworks: the single degenerate and the double degenerate scenarios. The two competing schemes for the 
galaxy formation, namely the monolithic collapse and the hierarchical clustering, are also taken into account, by 
considering the histories of star formation increasing and decreasing with redshift, respectively. We calculate 
j^ \ the Type la supernova rates through an analytical formulation which rests upon the definition of the SN la rate 

following an instantaneous burst of star formation as a function of the time elapsed from the birth of the progenitor 
system to its explosion as a Type la supernova (i.e. the delay time). What emerges from this work is that: i) we 
^ ■ confirm the result of previous papers that it is not easy to select the best delay time distribution scenario from the 

t:j- I observational data and this is because the cosmic star formation rate dominates over the distribution function of 

^O ■ the delay times; ii) the monolithic collapse scenario for galaxy formation predicts an increasing trend of the SN 

la rate at high redshifts (mainly due to the contribution by massive spheroids), whereas the predicted rate in the 
framework of a decreasing cosmic star formation rate, more in agreement with the hierarchical scenario, drops 
r^ I dramatically at high redshift; ///) for the elliptical galaxies we note that the predicted maximum of the Type la 

^r^ i supernova rate depends on the initial galactic mass. The maximum occurs earlier (at about 0.3 Gyr) in the most 

^— N I massive ellipticals, as a consequence of the assumed downsizing in star formation. In addition, we find that the 

Type la supernova rate per unit mass at the present time is higher in bluer ellipticals (i.e. the less massive ones). 
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1. Introduction 

The supernovae (SNe) of Type la are fundamental for 
understanding a number of astrophysical problems of 
primary importance, such as i) the SN progenitors, ii) 
the determination of cosmological parameters. Hi) the 
chemical enrichment of galaxies and iv) the thermal his- 
tory of the interstellar (ISM) and intracluster medium 
(ICM). The evolution of the rate of Type la SNe with 
cosmic time is a fundamental ingredient for the study 
of all these issues. The observed features of SNe la 
suggest that the majority of these objects may origi- 
nate from the thermonuclear explosion of a C-O white 
dwarf (WD) of mass ~ 1 .4 M© (Chandrasekhar mass) 
in binary systems (Chandra exploders). So each SN la 
should be the result of the explosion of the same mass. 
However, Phillips (1993) pointed out that there is a sig- 
nificant intrinsic dispersion in the absolute magnitudes 
at maximum light of local Type la SNe. This result 
was interpreted to arise from a possible range of masses 
of the progenitors or from variations of the explosion 
mechanism. 

Here we will consider only Chandra exploders for 
which two scenarios have been proposed: a) The Single 
Degenerate (SD) scenario, i.e. the accretion of matter 
via mass transfer from a non-degenerate companion, a 
red giant or a main sequence star (e.g. Whelan & Iben 
1973). In this scenario the mass range for the secondary 
componentsof the binary system is 0.8 -SM©, while the 
primary masses should be in the range 2 - SMq. The 
upper limit is given by the fact that stars with masses 
M > 8M0 ignite carbon in a non-degenerate core and 
do not end their lives as C-O WDs. The lower limit is 
instead obviously due to the fact that we are only in- 
terested in systems which can produce a Type la SN in 
a Hubble time. The clock to the explosion is given by 
the Ufetime of the secondary component, b) The Dou- 
ble Degenerate (DD) scenario, i.e. the merging of two 
C-O WDs which reach the Chandrasekhar mass and ex- 
plode by C-deflagration (e.g. Iben & Tutukov 1984). 
The merging is due to the loss of orbital angular mo- 
mentum due to gravitational wave radiation. In the Iben 
& Tutukov paper, the progenitor masses were defined 
in the range 5-8M0 to ensure two WDs of ~ O.TMq and 
then reach the Chandrasekhar mass. The clock for the 
explosion in this model is given by the lifetime of the 
secondary star plus the time necessary to merge the sys- 
tem due to gravitational wave radiation. This scenario 
requires the formation of two degenerate C-O WDs at 
an initial separation less than ~ 3/?o and this can occur 



by means of two different precursor systems: a close 
binary, and a wide binary. The two scenarios differ for 
the efficiency of the common envelope phase during the 
first mass transfer, and therefore for the separation at- 
tained at the end of the first common envelope phase. 

Different arguments can be found in favor or against 
both SD and DD scenarios and the issue of DD vs. SD 
is still debated (e.g. Branch et al. 1995, Napiwotzki et 
al. 2002, Belczynski et al. 2005). The Type la SN rate 
is the convolution of the distribution of the explosion 
times, usually called the time delay distribution func- 
tion (DTD), with the star formation history. Several at- 
tempts of comparing the Type la SN rate and the cosmic 
star formation rate evolution with redshift have already 
appeared (e.g. Gal- Yam & Maoz, 2004; Dahlen et al. 
2004; Cappellaro et al. 2005; Neill et al. 2006; Forster 
et al. 2006; Harris & Tonry, 2006; Poznansky et al. 
2007; Botticella et al. 2008, Blanc & Greggio, 2008; 
among others), but no clear conclusions arised yet on 
this point. 

In this paper, besides the above mentioned scenarios for 
Type la SN progenitors we will test the empirical DTD 
suggested by Mannucci et al. (2006), together with dif- 
ferent star formation histories. In particular, we will 
study the Type la SN rate in ellipticals and different cos- 
mic star formation rates, including a strongly increas- 
ing cosmic star formation with redshift, based on the 
monolithic scenario of galaxy formation and never con- 
sidered in previous studies. 

The paper is organized as follows: in Section 2 the 
adopted formulation for the SN la rate is presented, to- 
gether with the description of the adopted star formation 
histories. In section 3 the models for the distribution 
function of the delay times are presented and in section 
4 the predicted SN la rates are discussed and compared 
with the availale data. Finally, in section 5 some con- 
clusions are drawn. 

2. The computation of the Type la supernova rates 

From a theoretical point of view, the Type la SN rate 
is difficult to derive, because the nature of the progeni- 
tors of SNIa events is still an open question. 
One of the first model proposed for the calculation of 
the Type la SN rates was introduced by Greggio & Ren- 
zini (1983, hereafter GR83). This formulation is based 
on the Whelan and Iben (1973) model (SD scenario). 
Greggio (1996) revised the computation of the SN la 
rate in the framework of the SD model and suggested 
a more detailed criterion for the formation of a system 



which can eventually explode as Type la SN. Another 
model was proposed by Kobayashi et al. (1998) and 
Kobayashi, Tsujimoto & Nomoto (2000). They adopted 
the Hachisu et al. (1996; 1999) model and considered 
two possible progenitor systems: either a WD plus a 
red giant (RG) star or a WD plus a main sequence (MS) 
star, plus a metallicity effect delaying the formation of 
systems which could give rise to SNe la. Matteucci & 
Recchi (2001, hereafter MROl) calculated the Type la 
SN rate for different star formation histories in galax- 
ies in the framework of the SD scenario. They con- 
cluded that the best prescriptions to obtain SN la rates 
in agreement with the observations seems to be those 
of GR83 (with respect to the approach developed by 
Greggio 1996 and Kobayashi et al. 1998). Tornambe 
& Matteucci (1986) formulated a SN la rate in the DD 
scenario and applied it to galactic chemical evolution 
models. In the past years, many other authors (e.g. 
Ruiz-Lapuente & Canal, 1998; Madau, Delia Valle & 
Panagia, 1998; Dahlen & Fransson, 1999; Sadat et al., 
1998; Yungelson & Livio, 2000, Han & Podsiadlowski, 
2004; Scannapieco & Bildsten, 2005; Hachisu, Kato & 
Nomoto, 2008) presented formulations for the compu- 
tation of the rate of Type la SNe, alternative to that of 
GR83. Basically, most of them computed the SN la rate 
by introducing the function describing the delay time 
between the formation of the SN progenitors and its ex- 
plosion (i.e. DTD). One of the most recent formulation 
is that of Greggio (2005, hereafter G05) which will be 
adopted here. 

2.1. The formulation of the SN la rate ofGOS 

Following the G05 formalism, at a given epoch t the 
SN la rate is: 

Aja{t-T)l/,it-T)DTD(T)dT, (1) 



Xtr, 



where t is the delay time, defined in the range (t,, t^i), 
which are the minimum and maximum possible delay 
times ( i.e. the lifetimes of the secondaries in the SD 
model and the lifetimes of the secondaries plus the grav- 
itational time delays in the DD model), so that: 



I 



DTD(T)dT^ 1. 



(2) 



The constant ka is the number of stars per unit mass in 
one stellar generation and depends on the initial mass 
function (IMF): 

ka = I <t>(M)dM. (3) 



being (p{M) the chosen IMF, given by: 

(j){M) = CM-^^^^\ 
where C is the normalization constant derived from: 

• lOOMo 



(4) 
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M(f)(M)dM = 1. 



(5) 



In this work we adopt a Salpeter IMF, for which ka is 
equal to 2.83. Aja(t-T) is the fraction of binary systems 
which give rise to Type la SNe relative to the whole 
range of star masses (0.1 - lOOM©). In the formulation 
proposed by GR83 and MROl the realization probabil- 
ity Aja is instead relative to the mass range 3 - I6M0 
(i.e. the minimum and maximum total mass of the bi- 
nary system). In Eq. (1) A/a is a function of the epoch at 
which the stellar generation is born (f - t). However, in 
this work, this quantity is treated as a constant and it is 
chosen to reproduce the present day SN la rate observed 
value. This formulation shows that the basic ingredients 
for computing the SN la rates are the DTD and the star 
formation rate (SFR), tfr. 

For the analysis described in this paper we consider 
the SFR of typical elliptical galaxies and the cosmic 
SFR, namely the SFR in a unit comoving volume of 
the Universe containing galaxies of all morphological 
types. These functions are described in the following 
sections. 

2.2. The star formation rate in elliptical galaxies 

The SFR, i.e. the amount of gas turning into stars 
per unit time, is usually expressed in M© yr"'. In 
the case of elliptical galaxies we adopted the model 
of chemical evolution of Pipino & Matteucci (2004). 
The model is based on the assumption that ellipticals 
form by means of a rapid collapse of pristine gas, thus 
producing very high rates of SF (starburst-like regime). 
The SF is halted as the energy of the ISM, heated by 
stellar winds and SN explosions, exceeds the binding 
energy of the gas, thus producing a galactic wind. This 
happens on a timescale always less than 1 Gyr, vary- 
ing with the galactic mass. After the development of 
this wind, no star formation is assumed to take place. 
The SN feedback is taken into account together with 
the cooling of SN remnants. Massive but diffuse halos 
of dark matter around the galaxies are considered. The 
Salpeter IMF, considered constant in space and time, 
is adopted. The SFR has a simple form, given by the 
Schmidt (1959) law, with an efficiency of star formation 
which is higher in more massive objects which evolve 
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Figure 1 : The SFR per unit mass of typical elliptical galaxies of dif- 
ferent initial luminous masses, in units of Gyr"', taken from Pipino 
& Matteucci (2004). Magenta dot-dashed line: a model for IO'^Mq. 
Black solid line: a model for 10"Mo. Blue dashed line: a model 
for 10'~Mo. Note that the most massive ellipticals stop to form stars 
before the less massive ones. 



faster than less massive ones (inverse-wind scenario by 
Matteucci 1994, also called downsizing). The reason 
for this choice is that an increasing efficiency of star 
formation with mass reproduces very well the observed 
increase of [Mg/Fe] versus Mass in ellipticals. There- 
fore, the efficiency of star formation has been calibrated 
to reproduce such a relation (see Pipino & Matteucci 
2004 for details). As it can be seen in Fig. 1, we con- 
sider the specific SFRs (SSFR, i.e. the SFR per unit 
mass) corresponding to typical ellipticals of 10'", 10", 
and IO'^Mq of luminous mass. These histories of SF 
will then be used in the next chapter to predict the SN 
la rate under different assumptions about the DTD func- 
tion. 

2.3. The cosmic star formation rates 

The cosmic SFR is usually expressed in 
Moyr~'Mpc"^, in other words it is a SFR density 
(SFRd) in a unitary volume of the universe. Actually, 
the cosmic star formation history is measured up to z~6 
(see e.g. Hopkins, 2004; Hopkins & Beacom, 2006) 
thanks to the observations at UV, sub-mm and FIR 



wavelengths, although the error bars at high redshift 
are huge and the rates still uncertain. In the framework 
of a Lambda cold dark matter cosmology (ACDM) 
with Q.M = 0.3, Qa = 0.7 and Hq - 70km s^'Mpc^^ 
we have selected five different cases for the SFRd, as 
proposed in the literature. The models are those pro- 
posed by Calura Matteucci & Menci (2004, hereafter 
CMM04), whereas the others are from Madau, Delia 
Valle & Panagia (1998, hereafter MDP98), Strolger et 
al. (2004, hereafter S04) and Cole et al. (2001) and 
they are the best fits to various observations. It is worth 
noting that some of them predict a decrease of the SFRd 
at high redshift, in agreement with the hierarchical 
clustering scenario, but not supported by the present 
data indicating a rather constant behaviour (see fig. 
2). In particular, classical hierarchical clustering (HC) 
models for galaxy formation predict that ellipticals 
form continuously from the mergers of bulge-disk 
systems or other ellipticals, and that most galaxies 
never experience SFRs in excess of a few solar masses 
per year On the other hand, in the monolithic scenario 
the spheroids form first and the SFRs in massive 
systems can be as high as thousands solar masses per 
year CMM04 computed the theoretical cosmic SFRd 
as a function of redshift by means of detailed chemical 
evolution models for galaxies of different morphologi- 
cal types (Ellipticals, Spirals and Irregulars). Detailed 
descriptions of the chemical evolution models can be 
found in the works of Matteucci (1994) and Pipino 
& Matteucci (2004) for ellipticals (as described in 
sect. 2.2), Chiappini, Matteucci & Romano (2001) for 
spirals, and Bradamante, Matteucci, & D'Ercole (1998) 
for irregulars. 

The cosmic SFR density was then computed by 
CMM04 according to: 



M. 



Zivi 
PBj{z)(Y)B,ilpi(z), 



(6) 



where psj represents the B-band luminosity density 
(LD), (j;)bj is the B-band mass-to-light ratio and i/r, 
is the SFR for the galaxies of the /th morphological 
type. Therefore, in the CMM04 model, the total co- 
moving SFR density is given by the sum of the contri- 
butions of all the different morphological types ( see 
Calura & Matteucci, 2003 for details) and predicts a 
peak at the redshift of galaxy formation due to starbursts 
in spheroids (see Fig. 2; dot-dashed line). 

It is worth noting that the CMM04 cosmic rate is 
marginally consistent with the data for redshifts z > 3, 



and this is because it represents a real prediction, 
whereas the other rates are best fits to the data. On 
the other hand, the error bars of the data are very large 
especially at high redshift and it is not possible to 
distinguish whether the SFRd decreases, stays constant 
or increases for redshifts larger than 4. 

MDP98 presented two cosmic SFRd models (here- 
after MDPl and MDP2), expressed by two analytical 
fits, pmdp, and pmdpj: 



PMDP,(t) = adCe-'^'"^ + fl4(l - e-''/"^)] 



(7) 



and 



PMDP.iO = aie-'"'"" + fl4(l - e-"-'"') + astfe-'"'"' (8) 

in Mq yr"' Mpc"-'. The time tg is the Hubble time at 
redshift z, given by fg = 13(1 + z)~^^^K in Gyr and the 
values for the coefficients are given in Table 1. The 
MDPl model (Eq. 7, short-dashed line in Fig. 2) was 
built in order to fit the evolution of the observed comov- 
ing luminosity density at that time. The MDP2 model 
(Eq. 8, long-dashed line in Fig. 2) instead gives a larger 
SFRd at high redshifts and is designed to mimic the MC 
scenario (as in the model of CMM04). The last analyti- 
cal model we have considered is that computed by S04, 
assuming a modified version of the parametric form of 
the SFR as suggested by Madau et al. (1998), and it 
takes the correction for extinction into account (see Fig. 
2; dotted line). It is expressed in M© yr"^ Mpc"-' and is 
given by: 



psQ4(t) = ai(f^e '^"^ +046' 



d(t-to)/ay 



(9) 



where to - l3.47Gyr is the age of the Universe in Gyr 
(corresponding to z = 0), and the values of the other 
parameters are given in Tab. 1. In Fig. 2 all the de- 
scribed models are plotted together as a function of red- 
shift assuming two different epochs of galaxy formation 
Zf - 6, 10. An epoch of galaxy formation as large as 
Zf = 10 is justified by the fact that a number of well 
formed massive galaxies has been observed at z > 5 
(e.g. Mobasher et al. 2005; Vanzella et al. 2008). The 
conversion f - z is made by adopting the ACDM cos- 
mology, as defined before. The models are compared 
with the data provided by Hopkins (2004) and the para- 
metric form introduced by Cole et al. (2001) and used 
by Hopkins & Beacom (2006) to fit their SFR measure- 
ments (Fig. 2, solid line): Pi,(z) - (a + bz)h/(l +{zlcY), 
here with h=0.7. For a Salpeter IMF, the fit to the data 
gives a = 0.01334, /7 = 0.175, c = 2.93, J = 3.01 (see 
Blanc & Greggio, 2008). 



3. The distribution functions of the delay times 
(DTDs) 

For this work we have selected the formulation of 
MROl for the SD scenario (see also Matteucci et al. 
2006, hereafter M06), that of G05 for the DD wide 
and close scenarios and the empirical DTD proposed 
by Mannucci Delia Valle & Panagia (2006, hereafter 
MVP06). Different explosion models predict different 
delay times between the formation of the progenitor 
system and the SN explosion. In particular, the clock 
for the explosion in the case of the SD model (GR83 
and MROl) is given by the lifetime of the secondary 
star, ranging between ~0. 03-0. 04 Gyr (i.e. the lifetime 
of a 8Mq star) and 10 Gyr (i.e. the Hfetime of a O.8M0 
star This ensures that this model is able to predict a 
different from zero present time SN la rate for those 
galaxies where the star formation must have stopped 
several Gyr ago, such as ellipticals. In the DD sce- 
nario, the delay time is given by the sum of the life- 
time of the secondary plus the gravitational time delay 
(i.e. the time needed by the two stars to merge after 
loss of angular momentum due to gravitational wave 
emission, see Landau & Lifschitz 1962). In the DD 
scheme, the minimum delay time is given by the sum 
of the minimum nuclear delay time (0.03-0.04 Gyr) and 
the minimum gravitational delay time. This minimum 
gravitational timescale was suggested to be about 1 Myr 
by G05. Mannucci et al. (2005, 2006) suggested, on 
the basis of observational arguments, that there are two 
populations of the SN la progenitors, where a percent- 
age from 35% to 50% of the total SNe la explode soon 
after their stellar birth, namely inside ~ 0. 1 Gyr (the so- 
called prompt SNe la), while the rest has a much wider 
distribution of lifetimes, exploding over a long period 
of time (the so-called tardy SNe la). M06 tested this 
hypothesis in models of chemical evolution of galaxies 
of different morphological type (ellipticals, spirals and 
irregulars). They showed that this proposed scenario is 
compatible also with the main chemical properties of 
galaxies as long as the fraction of prompt SNe la is no 
more than 35%. Recently, a direct measurement of the 
DTD function has been reported by Totani et al. (2008), 
on the basis of the faint variable objects detected in 
the Subaru/XMM-Newton Deep Survey (SXDS). They 
concluded that the DTD function is inversely propor- 
tional to the delay time, i.e. the DTD can be well de- 
scribed by a featureless power law (DTDoc t"", with 
n ^1) in a range t =(0.1-10) Gyr. Such a DTD seems 
to support the DD scenario because the SD scenario, in 
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Figure 2: Cosmic SFR densities as functions of redshift, being Zf = 6 (upper panel) and zf = 10 (lower panel) the galaxy formation epoch. 
The SFR densities are in units of Mq yr"' Mpc"^. The models are from: Calura, Matteucci & Menci (2004) (dashed-dotted line); Madau, 
Delia Valle & Panagia (1998) (short-dashed line); Madau, Delia Valle & Panagia (1998) (long-dashed Une); Strolger et al. (2004) (dotted 
line). The different behaviours of the cosmic SFR are clearly visible toward higher redshifts. The data are from the compilation of Hopkins 
(2004) and the solid line is the best fit function to the data in the figure by Cole et al. (2001). This latter curve does not depend on the assumed 
galaxy formation epoch. 
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Table 1: Values of the coefficients for the analytical SFR densities. In the first column are indicated the models from Madau Delia Valle & 
Panagia (1998), Pmdpi ^^'ipMDPl, and the one from Strolger et al. (2004), pso4- 



some binary synthesis codes, predicts prominent peaks 
in the DTD at characteristic timescales (see e.g. Ruiz- 
Lapuente & Canal, 1998; Yungelson & Livio, 2000; 
Belczynski et al. 2005; Meng et al. 2008). However, 
other SD models predict other DTD shapes similar to 
the observed DTD (e.g. G05; M06; Kobayashi, Tsuji- 
moto & Nomoto, 2000). 

3.1. The DTD of MVP06 

MVP06 had used three important observational re- 
sults to derive, on empirical grounds, the DTD between 
the formation of the progenitor star and its explosion 
as a SN. The fraction of prompt SN la is suggested to 
be about 50%. This suggestion was based on the the 
enhancement of the SN la rate in radio-loud early-type 
galaxies, i.e. the fact that early-type radio-loud galaxies 
show a strong enhancement of the SN la rate with re- 
spect to the radio-quiet ones (see also Delia Valle et al. 
2005). 

The result is a bimodal DTD given by the sum of 
two distinct functions: a prompt Gaussian centered at 
5 ■ lO^yr and a much slower function, either another 
Gaussian or an exponentially declining function. An 
analytical formulation for this DTD is given in the 
workofM06: 



of the DTD of MVP06, a slope y = 2.0 should be as- 
sumed for systems with mass in the range 5 - 8M0 (i.e. 
in this mass range system where the primary and sec- 
ondary mass are almost equal are preferred), whereas a 
negative slope, y ^ -0.8, -0.9, should be adopted for 
masses lower than 5Mq (that means that systems where 
Ml » M2, in this mass range, are favored). 

Following M06, we assume, for this empirical DTD, 
the prescription of the SD model, thus t, and Tx are the 
lifetimes of a 8M0 and O.8M0, respectively. 



3.2. TheDTDofMROl 

The DTD of MROl is computed adopting the follow- 
ing formalism: 



DTD{t) oc 4>{M2)M2 



(11) 



that corresponds to the the SN la rate for an instanta- 
neous starburst. The function (j>{M2) is the mass func- 
tion of the secondary component, and in this case is 
given by: 



0(M2) = 2'+''(l+r)M^(M; 



-y tuA^^^^y") 



logDTD(T) 



1.4 -50(logT- 7.7)2 
-0.8 -0.9(logT- 8.7)2 



forT< 10^ 
forT> 10^ 
(10) 
where the time is expressed in years. This new formu- 
lation, being analytical, is easy to implement in a galac- 
tic chemical evolution code, in particular adopting the 
formalism developed by G05 (Eq. 1). As pointed out 
by M06, in the SD scenario, such a DTD can be justi- 
fied if one assumes that the function describing the dis- 
tribution of the mass ratios inside the binary systems, 
fifi) = 2^+^(1 + y)fi'>' (Tutukov & Yungelson, 1980), is 
a multi-slope function. In this expression, jj. = M2/MB 
is the mass fraction of the secondary and f{fi) is nor- 
malized to 1 between yU = and fi - 1/2. In the case 



■M^g '-'''>) /(-s-y), 
(12) 
93with s - 1 + X (x is the Salpeter index). The derivative 
93M2 = dM2/dt was obtained adopting the inverse of the 
formula (MROl): 

t(M) = 1 0*^' •^^'^^ Vl.79-0.2232(7.764-log M))/0. 1116 ^^3-) 

which defines the relation between the stellar mass (M) 
and the main sequence lifetime (t), with M expressed 
in M0 and r in yr 

Here the parameters are: 7 = 2 for a one-slope func- 
tion /(//), s - 2.35, Mh - max(2M2, Mmin), Mb - 
M2 + 0.5(M„«.,) with M,ni„ = 3Mo and M„a^ = I6M0 
i.e. the minimum and maximum masses of the system 
(see MROl, GR83). 



3.3. The DTD ofGOSfor the DD scenario 



4. Results 



In order to compute the DTD of the DD model, for 
both wide and close channels, G05 refers to systems 
with 2Mo < Mi,M2 < 8M0, from which most dou- 
ble C-O WDs form. Typically, the WD mass of both 
components ranges between 0.6 and 1.2 M©, so that the 
lifetime, t„, ranges between about 1 and 0.03 Gyr. The 
DTD depends on the distributions of both t„ and Tg„, 
with early explosions provided by systems with short t„ 
and Tg„. The delay time is given by t = t„ + Tg^^., with 
Tg^^. ranging between IMyr and a maximum value that is 
larger than t - t„ for all t„ (at least for a total delay time 
up to the Hubble time) in the wide DD scheme, whereas 
in the close one it is correlated with t„ The DTDs pro- 
posed by G05 for the wide and close DD schemes es- 
sentially depend on the following parameters (we refer 
the reader to G05 for details): 

• the distribution function of the final separations, 
which follows a power law: n(A) oc A'""; 

• the distribution of the gravitational delay times. It 
follows a power law, n(Tgv^.) oc t^^,; 

• the maximum nuclear delay time, t„^x- G05 treats 
Tn,x as a parameter, comparing the results obtained 
adopting 0.4, 0.6, 1 Gyr corresponding to the life- 
times of a 3, 2.5 and 2 M©, respectively. 

Following the results obtained by G05, we assume 



0.4Gyr, /3« = -0.9 and Pg 



-0.75. Then the 



DTD is given by: 



DTD(t) oc 



DTD(t) oc 



/-*mm 



•minijn.xj) 



n{T„)S (t, T„)dTn 



wideDD 

(14) 



■min{Tn.x,T) 



n{Tn)S (t, Tn)dTn cloScDD 

(15) 

complemented with DTD{t) - for t < t, and t > Tj. 
where S "^ and S '^ are functions of Pa and /3g, respec- 
tively (see G05 for details). We remark that, in the dif- 
ferent models, the number of prompt Type la SNe has 
an important impact in defining the shape of the DTD. 
In the DTD models proposed by MVP06, MROl (SD 
scenario) and G05 (DD scenario) about the 50%, 13% 
and 7% of the SNe la, respectively, are found to be 
prompt. 



By convolving the different DTDs with the various 
SFRs, described in the previous sections, we obtain the 
following results. We first concentrate on a particu- 
lar morphological type, ellipticals and then expand the 
analysis to the cosmic rates. 

4.1. The predicted Type la rates in ellipticals 

We consider elliptical galaxies of different baryonic 
mass (10'", 10", 10'^ Mq) adopting the SFRs of Pipino 
& Matteucci (2004). The parameter A/a, i.e. the frac- 
tion of systems which are able to originate a SN la ex- 
plosion, has been chosen a posteriori to ensure that the 
predicted present day SNe la rate in ellipticals is com- 
patible with the observed average value, given by Cap- 
pellaro et al. (1999): 0.18 + 0.06 SNu, being lSNu= 
l^A^/lO'^LoB/cenfMry. In particular, we have chosen 
the realization fraction A/a for each DTD, by reproduc- 
ing the present time Type la SN rate in a typical ellipti- 
cal of present time luminous mass of about 3.5- 1O'"M0, 
corresponding to an initial luminous mass of lO^'M© . 
To obtain the present time Type la SN rate, in units of 
SNe century^\ for an elliptical with a stellar mass of 
3.5- 10'"Mo, we multiplyied the Cappellaro et al. (1999) 
rate by the blue luminosity predicted by our photomet- 
ric model for such a galaxy (Lg = 4 ■ IO^Lb©), thus ob- 
taining a Type la SN rate of 0.072 SNe century^^ (see 
Table 2). 

Then the Aja obtained for this typical elliptical and 
a given DTD has been applied to ellipticals of different 
mass. In particular. Figs. 3, 4 and 5 show the results for 
an elliptical of initial mass IO'^Mq , IO^'Mq , IO'^Mq , 
respectively, versus redshift. We assumed z/ = 6 for 
the galaxy formation epoch. It is interesting to note that 
the maximum in the Type la SN rate depends on the 
galactic mass through the star formation history. In par- 
ticular, for a galaxy with initial mass of 10 '^M©, which 
suffers a very strong burst of SF with high efficiency and 
for a shorter time than less massive galaxies (downsiz- 
ing in star formation), the maximum of the Type la SNe 
occurs at about 0.3 Gyr, whereas in an elliptical with 
initial mass of lO^'M© the peak is at about 1 Gyr. This 
particular behaviour of the Type la SNR with galactic 
mass is fundamental to reproduce the increase of the 
average [Mg/Fe] in ellipticals, with galactic mass (see 
Matteucci 1994). It is worth noting that it is very diffi- 
cult to reproduce such a trend in the classical HC frame- 
work (see e.g. Thomas et al. 2002). 
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Figure 4: The predicted SN la rates for an elliptical galaxy of initial 
luminous mass of 1O"M0 versus redshift. The line types and units 
are the same as in Fig. 3. 



Figure 3: The predicted SN la rates for an elliptical galaxy of initial 
luminous mass of lO'^M© versus redshift . The results for differ- 
ent DTD models are shown: MROl (solid lines) and MVP06 (short- 
dashed lines) for the SD scenario; G05 for wide (dashed lines) and 
close (dotted lines) DD scenario. The Type la SNR is in unit of 
century"' and is shown on a logarithmic scale. All the predictions 
are normalized in order to reproduce the observed present day Type 
la SN rate in ellipticals (see text). Here, the redshift of galaxy forma- 
tion is assumed to be zf = 6. 







' 


, 1 1 1 . 1 1 1 . 


i 


10= 
10 






/€#-| 


-E 


1 


r 


/ 


'^ \ 


n 


0.1 




1 




1 



Figure 5: The predicted SN la rate for an ellitical galaxy of initial 
luminous mass of 10'~Mo, versus redshift. The line types and units 
are the same as in Fig. 3. 



In Table 2 are indicated the total number (SNIa) of 
SNe la explosions over the galactic lifetime, the pre- 
dicted SN rates (SNR) at the present time, in units of 
SNe century'^^ , and the ejected masses of Fe (after the 
galactic wind) for typical elliptical galaxies of initial lu- 
minous masses Mi of 10'", 10", lO'^M© , respectively. 
The total Fe mass produced by each galaxy was com- 
puted assuming that a single SN la typically produces 
about 0.6Mq of Fe (e.g. Nomoto et al. 1997). For the 
three different elUpticals, A/^ - 0.0019 , for the MROl 
model, A,„ = 0.012, for MVP06, A,a = 0.00021 for 
close DD and Aja - 0.0005 for wide DD are assumed. 
The computation of the Fe masses is useful in order to 
study the enrichment of the intracluster medium. The 
predicted values of the Fe masses are in agreement 
with those predicted by Matteucci & Vettolani (1988) 
and by Matteucci & Gibson (1995). We can note, in 
all the three cases, that the G05 DD models predict 
less SNe la and less Fe with respect to the SD mod- 
els. The last analysis we have made on ellipticals is 
based on the study of the SN la rate per unit mass, 
i.e. expressed in SNuM (Mannucci et al. 2005), be- 
ing ISNuM - IS N/century/lO^^^MQ. As suggested by 
Mannucci et al. (2005), the SN rate normalized to the 
stellar mass of the parent galaxies contains unique in- 
formation on the IMF of stars in the range of masses 
between about 3 and 100 Mq. Therefore, it is a very 
powerful tool for understanding the formation and the 
chemical evolution of galaxies and constraining their 
star formation histories. Here we focus the attention 
on the dependence of the SN la rate of elliptical galax- 
ies on their (B - K) color, which is an indicator of the 
mean age of the population. The SFR is related both 
to the morphology and to the color of the galaxies (e.g. 
Kennicutt 1998), but the SFR-color relation is proba- 
bly tighter than the SFR-morphology one as the color is 
more directly related to stellar populations than to mor- 
phology. Assuming that the stellar mass at the present 
time of the elliptical galaxies is about ~ 40% of the ini- 
tial luminous mass (see Pipino & Matteucci 2004), the 
predicted SNRs at the present time, in SNuM (see Fig- 
ure 6), for the three different galaxies are correspond- 
ing to the DTD models from MROl, MVP06 and G05 
wide and close DD. The colors have been computed by 
means of the spectrophotometric code of Jimenez et al. 
(2004). We can see, in Figure 6, that the bluer is the 
{B - K) colour of the parent galaxy, the higher is the 
predicted rate per unit mass at the present time (the first 
point on the left corresponding to an elUptical of initial 
luminous mass of lO'^M©). This behaviour can be ex- 
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Figure 6: The predicted Type la SN rate at tlie present time vs. (B-K) 
color of elliptical galaxies. The symbols are the predictions obtained 
by adopting the DTD from G05 in the wide scheme (stars), the DTD 
from G05 in the close scheme (open circles), the DTD from MROl 
(squares) and the DTD from MVP06 (triangles). The present day SN 
la rate is expressed in SNuM. i.e. in SNe per unit mass per century 
and it has been computed assuming that the present time stellar mass 
of the galaxy is about the 40% of its initial stellar mass. For a typical 
elliptical of initial luminous mass of lO" Mq we find a SNIa rate of 
about 0.018 SNuM. The (B - K) color increases with the mass of the 
galaxy, being reddest for a lO'^Mo-The very small dift'erence between 
the present time predicted values of the models for the galaxies with 
initial masses lO'" and IO'-Mq is due to the choice of the realization 
probability A/o, calibrated to reproduce the present time Type la SN 
rate of a 10"Mn. 



plained by the fact that the peak in the SFR (and hence 
the maximum in the SN la rate) occurs later in the less 
massive (bluer) ellipticals. 

Unfortunately, the available data give only an aver- 
age rate of Type la SNe as a function of morphological 
type and not as a function of colors for ellipticals of dif- 
ferent masses (for example SNuM = (0.027+°o[') for 



{B - K) > A.l in Mannucci et al. 
more recent Mannucci et al. 2008). 



2005, see also the 



4.2. The predicted cosmic SN la rate densities 

We present here the cosmic SN la rate densities, 
namely the number of Type la SN explosions per yr 
and per Mpc-', at low and high redshifts. The rates are 
computed by adopting the SFR densities described be- 
fore (Fig. 2). Here we compare the predicted cosmic 
Type la SN rates in order to define the differences due 
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Model 


SNla 


SNR(century-') 


Mf,(Mo) 




Ml 


= 10 '"Mo 




MROl 


4.33- 10' 


0.0108 


2.6-10' 


MVP06 


3.32- 10' 


0.0109 


1.99-10' 


G05 wide 


1.1-10' 


0.0107 


0.66 - 10' 


G05 close 


0.46 ■ 10' 


0.0105 


0.28 - 10' 




Ml 


= 10"Mo 




MROl 


2.7 ■ lO'^ 


0.072 


1.6- 10** 


MVP06 


2.23 ■ lO'^ 


0.072 


1.34 -10** 


G05 wide 


0.7 ■ 10^ 


0.072 


0.42 - 10** 


G05 close 


0.3 ■ lO"* 


0.072 


0.18 -10** 




Ml 


= lO'^Mo 




MROl 


1.86- W 


0.437 


1.12- lO** 


MVP06 


1.47- lO'' 


0.434 


0.88 - W 


G05 wide 


0.47 ■ lO'^ 


0.412 


0.28 - lO'^ 


G05 close 


0.2 ■ 10^ 


0.442 


0.12 -lO'' 



Table 2: The predictions for elliptical galaxies of initial luminous mass of 10 , 10 , 10 M©. In the first column the various DTD models 
are specified, while the second column shows the the total number of SNe la exploded over an Hubble time, predicted by each model. The 
predicted current SN la rates (SNR), in units of century"' , and the total Fe masses produced, in units of solar masses, are indicated in the third 
and fourth columns, respectively (see text). 



to the choice of different scenarios for the SN la progen- 
itors (namely, the SD and DD channels) and of differ- 
ent galaxy formation schemes (namely, the increasing, 
decreasing and constant SFRd). In all cases the realiza- 
tion probability A/^ was chosen a posteriori in order to 
reproduce the observed rate of Mannucci et al. (2005, 
see Table 3) at redshift z^0.03. Figs. 7, 8, 9 and 10 
illustrate the predicted cosmic Type la rates, separately 
for each DTD, compared with the observed ones at dif- 
ferent redshifts, assuming Z/ = 6 for the galaxy forma- 
tion epoch. The values of the parameter Aja are labeled 
in each panel, with line types refering to the different 
SFRds. Some observed Type la rates at different red- 
shifts are collected in Table 3. 

As it can be seen in these Figures, the predicted rates 
at low redshift (z < 1 ), are similar for all the differ- 
ent DTDs. At higher redshifts, the rates computed with 
the SFRds of MDPl (short-dashed Unes) and Cole et al. 
(2001, solid lines) show a steeper decline than the oth- 
ers, whereas those computed adopting MDP2 (i.e. for 
a MC scenario; (long-dashed lines) display a mild rise 
up to z ~ 4 . The most marked differences are shown 
when the theoretical SFRd of CMM04 is used. In fact, 
in the scenario of CMM04 the elliptical galaxies form at 
very high redshifts and suffer a quite strong star forma- 
tion, therefore producing large SN la rates at high red- 



shifts. The spirals instead contribute to the SFRd later 
and reach a peak at z ~ 2. This is a big difference with 
the models which follow a hierarchical scheme, in other 
words with models predicting a decrease of the cosmic 
SFRd at high redshifts and assume that ellipticals form 
mainly at lower redshift (z ~ 1 - 2). As a consequence 
of this, the detection of Type la SN rates at high redshift 
(z > 3) would help to decide which scenario of galaxy 
formation is favoured. In particular, the models adopt- 
ing the SFR computed from CMM04 for a MC scenario 
are the only ones predicting a very high SN la rate at 
high redshift. The curve obtained by convolving the 
bimodal DTD of MVP06 and the SFR of CMM04, as- 
suming Zf - 6, shows a high peak at early times (high 
redshifts), predicting a larger number of Type la SNe 
(~ 73.3 - 10-"* SNe la yr-lMpc-^ at z ~ 4.8) with re- 
spect to the other models. It is important to point out 
that, at high redshifts (z > 1), the majority of the models 
seem to overestimate the data by Dahlen et al. (2004) 
solid squares). It has been suggested that these data are 
very difficult to reproduce by means of standard SN rate 
models (Dahlen et al. 2004, Mannucci et al. 2005), un- 
less a large delay time (~ 4 Gyr) between the epoch 
of star formation and the explosion of Type la SN is 
assumed. This means that one should wait a time com- 
parable to ~ 4 Gyr after the beginning of star formation 
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z 


SN rate 
(10 *Mpc 3yr 1) 


Reference 


z 


SN rate 
(10 *Mpc 3yr ') 


Reference 


0.03 


0.28±0.11 


Mannucci et al. 2005 


0.55 


0.52±0.14 


Pain et al. 2002 


0.10 


0.32± 0.15 


Magdwick et al. 2003 


0.55± 0.05 


2.04± 0.38 


Barris & Tonry 2006 


0.11 


0.37± 0.10 


Strolger et al. 2003 


0.65± 0.05 


1.49±0.31 


Barris & Tonry 2006 


0.13 


0.20± 0.08 


Blain et al. 2004 


0.75± 0.05 


1.78±0.34 


Barris & Tonry 2006 


0.2 


0.1 89± 0.042 


Horesh et al. 2008 


0.75± 0.25 


0.43±+0.36-0.32 


Poznanski et al. 2007 


0.24± 0.05 


0.17±0.17 


Barris & Tonry 2006 


0.8 


1.57±+0.5-0.4 


Dahlen et al. 2004 


0.3 


0.34±+0.16-0.15 


Botticella et al. 2008 


1.2 


1.15±+0.7-0.5 


Dahlen et al. 2004 


0.35± 0.05 


0.53± 0.24 


Barris & Tonry 2006 


1.2+0.2 


0.75±+0.35-0.30 


Kuznetsova et al. 2008 


0.4 


0.69±+0.2 - 0.4 


Dahlen et al. 2004 


1.25±0.25 


1.05±+0.45-0.56 


Poznanski et al. 2007 


0.4± 0.2 


0.53±+0.39-0.17 


Kuznetsova et al. 2008 


1.55±0.15 


1.2±+0.58-0.12 


Kuznetsova et al. 2008 


0.45± 0.05 


0.73± 0.24 


Barris & Tonry 2006 


1.6 


0.44±+0.5 - 0.3 


Dahlen et al. 2004 


0.46 


0.48± 0.17 


Tonry et al. 2003 


1.75±0.25 


0.81±+0.79-0.60 


Poznanski et al. 2007 


0.47 


0.42± 0.06 


Neil et al. 2006 









Table 3: A compilation of the available observational cosmic Type la rates at different redshifts up to z~1.75. The SN la rates are expressed 
in units of lO^^Mpc^^yr"' . These data are compared with the predicted cosmic rates in Figs. 7-15 



in order to observe a significant contribution from Type 
la SNe to the chemical enrichment of any astrophysi- 
cal system. However, chemical evolution models indi- 
cate that in the solar neighborhood the time at which 
the Fe production from SNe la starts to become im- 
portant is ~ IGyr after the beginning of star formation 
(Matteucci & Greggio 1986; MROl), well reproducing 
the [a/Fe] values as a function of the [Fe/H] ratio ob- 
served in Galactic field stars. Moreover, MVP06 have 
shown that a large fraction of SNe la should arise from 
systems exploding on timescales on the order of 40Myr, 
to explain the Type la SN rates observed in radio-loud 
elliptical galaxies. The data observed by Dahlen et al. 
(2004) were the first collected at redshifts > 1 for the 
Type la SN rate and are likely to represent lower limits 
to the actual values (and are still uncertain). However, 
recently Kuznetsova et al. (2008) presented a new mea- 
surement of the Type la rate up to redshift of 1 .7, using 
two samples collected by the Hubble Space Telescope. 
They applied a novel technique for identifying Type la 
SNe, based on a Bayesian probability approach, and 
their results (open circles) appear to be slightly larger 
than those of Dahlen et al. (2004). In this way, the 
agreement between the predictions of the models pre- 
sented here and the observations improves and there is 
no more need to assume very large delay times. More- 
over, it is very important to note that the assumed star 
formation history dominates over the assumed DTD and 
that one can obtain a cosmic SNIa rate bending down at 
z~2 with a drcreasing SFRd and any DTD. Therefore, 



the suggestion of the long delay time for the SNIa rate 
is not a robust conclusion anyway. 

In Figs. 11, 12, 13 and 14 we show the results ob- 
tained by assuming a different epoch of galaxy forma- 
tion, Zf - 10. The values of the realization probability 
Ala are the same assumed in the case with z/ = 6 (as 
labeled in the figures). We note in these figures that the 
models with Zf = 10 , computed in the framework of 
a decreasing SFRd with redshift, show a SNela peak at 
a slightly larger z compared to the models with z/ = 6, 
and that all the models show larger SNIa rates at z ~ 
2. The behaviour of the models in the redshift interval 
< z < 1 (the range better covered by observations) 
does not vary appreciably when changing the redshift 
of galaxy formation. Finally, in Fig. 15 we show the 
results obtained with the cosmic SFRd of Cole et al. 
(2001) coupled with all the explored DTDs. Here the 
agreement looks good for almost all the DTDs, with 
perhaps the exception of the DTD relative to the close 
DD model, thus confirming the conclusion that it is very 
difficult to select a DTD from the observed cosmic Type 
la SN rate. 



5. Conclusions 

We calculated the cosmic SN la rate densities (i.e. 
the rate per unit comoving volume) and the rate of the 
explosion of SNe la in typical elliptical galaxies adopt- 
ing the formalism proposed by G05. This formulation 
rests upon the definition of the SN la rate following an 
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Figure 8: As in Fig. 7, but for tlie MROl DTD. 



Figure 7: The cosmic SN la rate density for the MVP06 DTD. This 
DTD is convolved with the SFRd from CMM04 (dot-dashed line), 
MDP2 (long-dashed line), MDPl (short-dashed line), S04 (dotted 
line) and Cole et al. (2001, solid line). The short-dashed, dotted 
and solid curves show a decrease at redshift z > 1 with respect to 
the other models (dot-dashed and long-dashed lines). The SN la rate 
is expressed in units of 10"''yr"'Mpc"''. The redshift ranges from 
to 6 (i.e. the epoch at which the star formation is supposed to have 
started). The solid squares are the data by Dahlen et al. (2004), the 
open pentagon is by Botticella et al. (2008), the open circles are from 
Kuznetsova et al. (2008), the open triangles are from Poznanski et 
al. (2007), the open squares are from Barris & Tonry (2006), all the 
other data, given in Table 3, are represented as solid triangles. The 
predicted present day rates are normalized to reproduce the observed 
value at z=0 (see text). 
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Figure 9: As in Fig. 7, but for the wide DD DTD of G05. 
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Figure 10: As in Fig. 7, but for the close DD DTD of G05. 



Figure 12: The cosmic SN la rate density for the MROl DTD assum- 
ing zf = 10 for the redshift of galaxy formation. The models, units 
and data are the same as in Fig. 7. 
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Figure 1 1 : The cosmic SN la rate density for the MVP06 DTD as- 
suming Zf = 10 for the redshift of galaxy formation. The models, 
units and data are the same as in Fig. 7. 



Figure 13: The cosmic SN la rate density for the wide DD DTD as- 
suming Zf = 10 for the redshift of galaxy formation. The models, 
units and data are the same as in Fig. 7. 
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Figure 14: The cosmic SN la rate density for the close DD DTD 
assuming zf = 10 for the redshift of galaxy formation. The models, 
units and data are the same as in Fig. 7. 
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instantaneous burst of star formation as a function of 
the time elapsed from the birth of the progenitor system 
to its explosion as a Type la SN (i.e. the delay time). 
This function has been termed as DTD and accounts 
for the SNe la progenitor scenarios and for the initial 
mass function (IMF). Different DTDs and histories of 
SF have been considered. In all cases, a Salpeter IMF 
was adopted. For all galaxies it was assumed a galaxy 
formation epoch at redshift zy = 6. The reason for chos- 
ing Zf = 6 is given by the fact that the SFR density has 
been measured up to this redshift (Cole et al. 2001). We 
have also shown the results obtained assuming z/ = 10 
as epoch of galaxy formation. The evolution of the cos- 
mic SN rate with redshift contains, in principle, unique 
information on the star formation history of the Uni- 
verse, the IMF of stars, and the Type la SN progen- 
itors. These are essential ingredients for understand- 
ing galaxy formation, cosmic chemical evolution, and 
the mechanisms which determined the efficiency of the 
conversion of gas into stars in galaxies at various epochs 
(e.g. Madau et al. 1996; Madau, Pozzetti, & Dickinson 
1998; Renzini 1997). We computed the cosmic Type la 
SN rates for different cosmic SFR histories. Our main 
result is the prediction of the expected number of explo- 



0.01 



Figure 15: The predicted cosmic SN la rates adopting the SFR density 
from Cole et al. (2001). The short-dashed line is for the MVP06 DTD, 
the solid line is for the MROl DTD while the dashed and dotted lines 
are for the wide DD and close DD DTDs, respectively (G05). The 
figure underlines the differences due to the choice of the DTDs. The 
data are the same as in the previous Figures. Note that this diagram is 
independent of the redshift of galaxy formation, since the SFR density 
adopted does not depend on zf- 
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sions of Type la SNe at high redshift (z > 2). 
Here we summarize our conclusions in detail: 

• we analyzed the effects of the various parame- 
ters entering the computation of the SN la rate 
and concluded that the realization probability Aja 
(the actual fraction of systems which are able to 
give rise to a Type la SN) should be in the range 
10"^ - 10"^, where the lower value is appropriate 
for the DD model. This quantity, in all models, is 
a free parameter and it was chosen by reproduc- 
ing the present time Type la SN rate in galaxies. 
Other parameters also play a key role in the com- 
putation of the Type la SN rate and these are: i) 
the mass range for the secondaries; ii) the mini- 
mum mass for the primaries; Hi) the efficiency of 
the accretion and iv) the distribution of the separa- 
tions at birth. Actually, all the different channels 
could contribute to the SN la events, each with its 
own probability (A/a). 

• We computed the Type la SN rate for specific ellip- 
tical galaxies of different initial luminous masses 
(10'"Mo, 10"Mo, IO'^Mq) for all the studied 
DTDs and for the SFRs suggested by chemical 
evolution models which best reproduce the char- 
acteristics of local ellipticals (Pipino & Matteucci 
2004). All the DTDs predict an early maximum in 
the Type la SN rate, between 0.3 and 1 Gyr, ac- 
cording to the galaxy mass. In particular, the max- 
imum predicted at about 0.3 Gyr corresponds to a 
typical elliptical of initial mass of lO'^M© while, 
for an elliptical of initial mass of lO'^Mg the peak 
is at about 1 Gyr. The dependence of the maxi- 
mum on the initial galactic mass is due to the ef- 
ficiency of star formation, which is assumed to be 
higher for the most massive ellipticals . We have 
also considered the predicted Type la SN rates per 
unit mass (SNuM) at the present time vs. (B - K) 
color relations for the three different ellipticals. 
We found that the bluer is the color (and hence the 
lower is the galactic mass), the higher are the pre- 
dicted SN la rates in SNuM at the present time. 
Unfortunately, a real comparison with data is not 
possible since the data for ellipticals of different 
masses are not available. 

• We computed the cosmic Type la SN rate in a uni- 
tary volume of Universe by adopting different cos- 
mic SFRds, as predicted by both models and ob- 
servations (Hopkins & Beacom, 2006). We found 



that the SF history largely dominates over the as- 
sumed DTD in the calculation of the Type la SN 
rate. Therefore, unless we know the cosmic star 
formation history, we cannot safely decide which 
DTD is better on the basis of the observed cos- 
mic SN la rate, in agreement with previous works 
(see also Forster et al. 2006; Blanc & Greggio 
2008, Botticella et al. 2008). In particular, it is 
not possible to predict the delay time for the explo- 
sion of SNe la on the basis of the cosmic star for- 
mation rate, since for example, in the hierarchical 
framework, it is predicted a decrease of the cosmic 
Type la SN rate for z > 1 irrespective of the cho- 
sen DTD, even for that of Mannucci et al. (2005; 
2006) DTD, where the fraction of prompt Type la 
SNe is 50%! 

The cosmic Type la SN rates for different SFRds 
(increasing, decreasing, constant with redshift) 
differ mostly for redshift z > 1 . We compared our 
results with the available data, although the highest 
redshift points are still very uncertain. Therefore, 
it was not possible to decide which scenario should 
be preferred, but high redshift data will greatly 
help to draw conclusions on this point. In particu- 
lar, we predict that the Type la SN rate for z > 2 
should be very high if the monolithic scheme for 
the formation of ellipticals is assumed and even 
higher if combined with the empirical DTD of 
Mannucci et al. (2005), as opposed to models with 
decreasing or constant SFRds, which all predict a 
significant drop of the Type la SN rate for any cho- 
sen DTD. Future observations with JWST of high 
redshift SNe will help in sheding light on this sub- 
ject. 
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